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1 . introduction 

Rainfall is a key link in .he hydropic cycle t «Kj I is • 

rrimarv heat source for the atmosphere. , 

distribution of latent-heat release, wh.cn 
accompanied by rainfall, modulates the large-scale 
circulations of’ihe tropics and in turn can impact 
midlatitude weather. This latent heat release is a 
v s-vf n }v'.cp changes between vapor, ntju.d. 

r h ", ” d 

'climate models can simulate cloud 

onlv crudely, thus reducing their confidence 

predictions on both global and regional scales. 

In this paper, NASA Tropical R ainfall Measuwn^ 
m; dop (TRMM) precipitation ladai t?R) 

“ XlUnform u on and d.c Goto, cl Coa.ec, ivc and 
Sfrc,™ Healing (CSH) algorithm are used to 
estimate the four-dimensional structure of g'oc.al 
monthly latent heating and rainfall profiles over the 
global tropics from December 1997 to October 200J 
Rainfall, latent heating and radar reflectivity str.icturt^ 
between ENSO (1997 -1998 winter) and non- .a IS 
( 1998-1999 winter) periods are examine an 
compared. The seasonal variation of heating over 
various geographic locations (i.e. Indian ocean vs west 
Pacific; Africa vs S. America) are also ana1y"d. In 
addition, the relationship between rainfaii, latent 
heating (maximum heating level), radar retlectivi y a 
SST are examined. 

2. Goddard Convective and Stratiform 
Heating (CSH) algorithm 

Diagnostic budget (House, ]9S2: Johnsoi 1 ,^1984) and 

cloud modeling studies (see a review by - 

have shown that the distribution of heating in the an il 
pm ion of tropical mesoscale cloud systems .s 

considerably different from the vertical profile of 

heatins in the convective region. Generally, 
evaporative cooling in the lower troposphere is 
dominant in the stratiform region in a.l ^rnilated 
convective svslems. On the other hand, 
condensation/deposition heating is dominant in the 
convective region in these convective systems. Has 
on these findings, a conveclive-siraiifoim h.at e 


aleorithm has been developed (Tao er at. 1 993). I hem 
are many sets of normalized heating profiles; each ^ 
has a convective profile and a stratifoim P' J 'T 
simulated by cloud resolving models or oiagnosv 
from budget calculations Each set represents di - a 
system oreanizalions as well as various geograpim- 
locai ions. Figure l shows the schematic diagram w 
the CSH algorithm. 
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Fio 1 Diagram showing the procedure Jot tie 
° latent heating profile from ^ 

heating (CSH) algorithm. The healing ptOji.e U 
normalized hr the surface precipitation rale F 0 >. 

Vertical profiles of latent beat release and their 
retrieval using the CSH algorithm associated ,n\» 
--veral TOGA CO ARE convective active pet lodsau 
examined by Tao ^(2000^ Tlte inputs 

(synthetic TRMM PR) derived rainfall and ^ufonn 
Z Li Results indicated that temporal variability ot 
CSH algorithm* retrieved latent healing profile* using 
radar estimated rainfall and stratiform amount n t 
cood agreement with that diagnost.ca y ^ 

However less rainfall and a smaller sttatifoiin 
percentage estimated by radar resulted m weaK,r 



(underestimated) latent heating profiles and lower 
maximum latent heating levels compared to those 
determined diagnostically. Rainfall information from 
SSM/I can not retrieve individual convective events 
due to poor temporal sampling. Time averaged heating 
profiles derived from SSM/L however, are not in bad 
agreement with those derived by soundings. 

Tao etat. (2001) used the TRMM rainfall information 
to estimate the four dimensional latent heating 
structure over the global tropics for one month 
(February 1998). The latent heating profiles derived 
from the CSH algorithm were compared to those 
derived from two other different heating algorithms, 
the Goddard Profiling (GPROF) heating, and the 
Hydrometeor heating (HH). The horizontal 
distribution or patterns of latent heat release from the 
three different heating retrieval methods are quite 
similar. They all can identify the areas of major 
convective activity [i.e., a well defined Intertropical 
Convergence Zone (ITCZ) m the Pacific, a distinct 
South Pacific Convergence Zone (SPCZ)] in the global 
tropics. The magnitude of their estimated latent 
heatine release is also in good agreement with each 
other and with those determined from diagnostic 
budget studies. However, the majoi difference among 
these three heating retrieval algorithms is the altitude 
of the maximum heating level. The CSH algorithm 
estimated heating prclilcs only show one maximum 
heating level, and the level varies between convective 
activity from various geographic locations. These 
features are in good agreement with diagnostic budget 
studies. 

3. TRMM PR Rainfall Data 

The TRMM Precipitation Radar (PR) operates at a 
frequency 13.8 GHz. This is a moderately attenuating 
frequency so that one of the major objectives of the PR 
algorithms is to estimate and correct for the path 
attenuation. Because the PR is a single-wavelength, 
single-polarization, non-Doppler radar, there arc only a 
few methods available to estimate rainfall. These 
include the Hitschfeld-Bordan (FIB) method, the 
surface reference technique or SRT, a ’’hybrid' 1 of 
these, and the mirror-image technique. The rain rate 
estimation algorithm uses a hybrid of the Hitschfeld- 
Bordan method and the SRT to correct for attenuation 
and derive an estimate of the range profile of the radar 
reflectivity factor, dBZ. The rain rate profile is ihen 
calculated from the dBZ profile and an appropriate Z-R 
relationship. The algorithm also includes surface 
clutter rejection and an attempt to correct for effects ot 
non-uniform beamfilling. 

Convective-stratiform classification of the rain is based 
on a combination of the vertical and horizontal 
structure of the radar reflectivity field. The vertical 


profile is checked for the presence of a brigHt-band 
(melting layer) by considering the behavior of the 
second derivative of the radar range profile. Unless the 
maximum dBZ exceeds a threshold, the presence of a 
well-defined bright-band is used to indicate stratiform 
rain. In cases where a clearly defined melting layer is 
absent, the horizontal rain structure is examined by 
means of a modified version of an algorithm designed 
for the analysis of ground-based radar data. 

Statistics of the instantaneous, high -resolution rain 
rates are compiled on a monthly basis over 0.5° x 0.5^ 
latitude-longitude grids. Near-surface rain rates and the 
rain rates at 3 altitude levels are stored according to 
rain type. The statistics include means, standard 
deviations and histograms of the rain rate, radar 
reflectivity factors, bright-band and storm heights.^ 
Please see Tao et cii (2001) for a brief description oi 
the PR rainfall algorithms and convective stratiform 
separation. 

4. Results 

Figure 2 shows three yearly mean latent heating at 
three different altitudes (2. 5 and 8 km) over the global 
tropics from the CSH algorithm using PR derived 
rainfall products and latent heating profiles in the iook- 
up table representing general tropical oceanic and land 
reeions (obtained by averaging profiles from the look- 
up table). The horizontal distribution pattern of the 
CSH estimated latent heating structures is very similar 
to the pattern of surface rainfall (not shown), especially 
at middle and upper levels. For example, a well 
defined ITCZ in the east and central Pacific and in the 
Atlantic Ocean, a distinct S. Pacific Convergence Zone 
(SPCZ). and broad areas of precipitation events spread 
over the continental regions, are present. Also, strong 
latent heat release (5 K/day or greater) in the middle 
and upper troposphere is always associated with 
heavier surface precipitation. Heating in the upper 
troposphere over the Pacific and Indian Oceans covers 
a much broader area than the healing over Africa, S. 
America and the Atlantic Ocean. The differential 
heatine distribution between land and ocean in the 
upper troposphere could generate strong horizontal 
gradients in the thermodynamic fields and interact with 
the global circulation. 

One interesting result from Fig. 2 is the relatively weak 
heating and cooling (-Mol K/day) at the 2 km level 
over the Pacific Ocean. This result-may be due to the 
fact that the moisture content is high over the Pacific 
ocean. Cooling by evaporation of raindrops in the 
lower troposphere could be weak over moist areas. 
Another possibility is that strong heating and cooling 
associated with the life cycle of many convective 
systems or seasons may cancel each other. 




fjo 2 Three year mean latent healing profiles at (a) 8, (b) 
5 and ( c ) 2 km (above ground level, AGL) over the 
»iobal tropics. Latent heating profiles representing 
tropical oceanic and general land regions from the 
CSH look-up table (LUT) were used. 

Time (monthly) series of latent heating profiles over 
the tropics and their deviations (from the mean) aie 
shown in Fig. 3 . The level of maximum healing is at 
7.5 km. The variation of the maximum heating level as 
well as the magnitude are quite small (Fig. 3a). 
However, there "arc cold anomalies during ENSO 
(1997-1998 winter) and warm anomalies during non- 
ENSO (1998-1999 winter) periods. These features are 
due to the fact that the PR observed a higher 
percentage of stratiform precipitation during the ENSO 
period which is conducive to retrieving stronger low 
level cooling compared to the non-ENSO. period. 
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